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Abstract 

We compute the survival probability of an initial state, with an energy 
in a certain window, by means of random matrix theory. We determine its 
probability distribution and show that is is universal, i.e. caracterised only 
by the symmetry class of the hamiltonian and independent of the initial 
state. 



In classical mechanics, temporal chaos is caracterised by the extreme sensi- 
bility of a trajectory to variation of initial conditions. No direct analog of this 
phenomenon has been found in quantum mechanics so far. On the other hand, 
numerical evidence has been accumulated [1], showing that energy levels of a 
quantum system, whose classical counterpart is chaotic, have a statistical be- 
havior described by Wigner's random matrix theory (RMT), on the mean level 
spacing scale. The question we want to adress is the following: are there specific 
predictions of RMT for quantum dynamics, which would caracterise the temporal 
behavior of " chaotic" quantum systems. 

We consider the following situation: The system is prepared in an initial state ip 
at time 0, with an energy in a certain window, centered at e and of width 2sl(e), 
where /(e) is the mean level spacing, and we want to compute the probability to 
find our system again in the state at a later time t. This quantity that we call 
the survival probability R is given by 



R 



(^|expzf# P{A)<p) 



(ip, P(A)ip) 

H is the hamiltonian of our system and -P(A) is the spectral projector on A. 
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We have chosen to take an energy in a range of the order of the mean level 
spacing in order to look at properties of the system which are independent of 
specific details. If (Xj, tpj) denote respectively the j l h eigenvalue and eigenvector 
of the Hamiltonian H, then the survival probability can be written as 



R 

where 



Ej=i VjXj exp 2-kitx, 2 



© Exi-M ( 2 ) 



Ej=i VjXj 

Vj = \(<P> ( 3 ) 



and if we define Xj by the relation 

Xj = e + Xjl(e) (4) 

{1 if x ■ ^ s 
otherwise ^ 

The Heaviside function ensures that there is at least one eigenvalue in A. 
What appears naturally in this expression is the time measured in units of the 
Heisenberg time 

so that 

If we look at this problem from the point of view of RMT, we will replace the 
Hamiltonian by a large N x N self-adjoint matrix, whose probability distribution 
is basis independent and therefore of the form 

e -w(X 1> ...,x N ) dH (g) 
Wigner's gaussian model corresponds to the choice 

AT N 

^=yEA? (9) 

Z 3=1 

The first conclusion to be drawn is that the survival probability is statistically 
independent of the initial state tp. This follows from the fact that the variables 
{yj}^ =l have a probability distribution, independent of ip and given by: 

1 / N \ N l 

My)dy = s\ E vj - 1 ] II v] d y (10) 



N 



\3=1 



The parameter (3 = 1, 2, 4 caracterise the symmetry class of the Hamilto- 
nian, respectively orthogonal, unitary and symplectic. Equation follows 
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easily from the Haar measure on the corresponding groups. Cat is a normalising 
constant. 

The variables {xj}J =1 are statistically independent of the variables {yj}^ =1 
and have a distribution given by 

-^exp-W (e + x/(e))A /3 (x)dx (11) 
D N 

where the Van der Monde determinant 

A(x)= J] (12) 

l<i<j<N 

comes from the change of variables Hij — > ipj)J =1 [2]. Dn is a constant of 
normalisation. 

We can take /(e) = ^ * , where p(e) is the density of states when N = oo. 
The problem that we need to solve now is to find the probability distribution of 
the survival probability p(K)dR in the N = oo limit. We find that R is not self- 
averaging i.e. p(R) is not a delta distribution concentrated on the mean value 
of R. On the other hand its probability distribution p(R) is universal, i.e. it 
depends only on the symmetry parameter (3, at least for a large class of W. 
There are two formulas for p(R), one more appropriate to small windows, another 
one to large windows. 
In the first case, we decompose p(R) into 

oo rp 

pW = Er>( fl W ( 13 ) 

n=l 1 - h 

where E n is the probability to find exactly n eigenvalues in A and p n {R\r) is the 
conditional probability density of R knowing that there are exactly n eigenvalues 
in A. 

It can be expressed as 



/s ^ poo 
E(x 1 ,...,x n ) d n x fin fa,..., z n ) d n zS 



( 



R- 



V 



Zj exp 2irirXj 

3=1 

(14) 

E(x u ... t x n ) = E{xu „--' Xn) (15) 

E n = J E(xi,...,x n )d n x (16) 
E (xi, . . . , x n ) being the probability density of finding the n eigenvalues in A at 

[X\ , . . . , Xn). 



2 \ 



J 
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Useful expressions for E (x±, . . . ,x n ) and E n can be found in [2] and [3]. It is 
expressible in terms of a determinant 



E (xi , . . . , x n j 



where 



det L/3 (xi\xj) 



E (l...n) 



(17) 



Kp is an operator whose kernel in the simplest case (3 — 2 is given by 

sin 7r(x — y) 



(x\y) 



ti(x - y) 



(19) 



defined on L 2 (— s, s). 

Universality comes from the fact that E (xi, . . . , x n ) is expressible in terms of 
the correlation functions and the latter ones depends only on (3, for a large class 
of W. W modifies only the density of states and therefore the mean level spacing 
/(e). 

This expression for p(R\r) is mostly useful in the small window limit, because 
when s —>■ 

E n ~ S* 1 



^ 2 +n(l-f) 



Moreover in this case we have 



lim s n E (sxi , . . . , sx n ) = A n IT 

s^O . . . . . 

i<i<j<n 



(20) 
(21) 



so that the probability distribution of R shows a scaling behavior 



lim s 
1 



-0-1 



pr 



1 - R 



> x 



TITS 



gp(X)d\ 



(22) 



the function ^(A) being given by 



(3-1 



gp{\) = Ap\ 2 



-Vl - A + - In A - In 1 + Vl - A 
2 2 



(23) 



On the other hand, one can see from eq (0) and (14) that the probability 
distribution of R is well defined at infinite times. Namely 



p n {R\r) =p n (R\oo) + 
as can be seen by an integration by parts where 



(24) 



Pn (R\00) 



i • ■ ■ i *n; 



2vr « 

n 



2vr 



/ 


n 






E 


] 


V 


j=i 





(25) 



Using an integral representation for the delta appearing in the definition ( |i0|) of 
the fi n , we can reexpress fl25|) as 



p n (R\oo) 



/ due u / drrJo(VRr) / dze~ uz J [ 

4:7lC n Je-ioo JO V ' [Jo 



rz) z 2 ' 



(26) 

e being any positive number, and Jq(x) the Bessel function. This expression can 
be simplified, considerably when /3 — 2, 4. 
In the unitary case ((3 = 2) one finds 



p n (R\oo) 



n — 1 



1 - R) 



71 — 3 

2 



(27) 



For a large window of energy, it is more appreciate to find another expression for 
p n (R\t). It is given as some integral over a Fredholm determinant Q. 
Q is a generating function for the variables {yj} and {xj} appearing in the defi- 
nition of R, eq (0). 



G{V]<p;\) = ^im/exp-)^^ f|X| [V cos (e7rT§| + tp) + $ 



(2* 



It can be expressed in terms of the operator Kp appearing in eq (|l^), when 

/3 = 1, 2 as 



Q = S, det ( oo + /C/j} e 



with 

E = [det (1-^)]^ 
and (7 is the multiplication operator by the function 
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1 + -!-[z + r cos (2tttx + </?)] 



(29) 
(30) 

(31) 



When the window is large (s >> 1) we can expand the determinant in powers 
of Kp, the first two terms of this expansion dominating the other ones [4]. 
One finds that the probability distribution is exponential. 



lim —p I — 



1 R 

exp 



a t 



a r 



(32) 



In the orthogonal case ((3 = 1), for example 
a{r) = 



4 — 2 |t| + |t| In 1 + 2 |t| 
2 



I I 1 2|t|+1 



if \r\ < 1 
if M > 1 



(33) 



One can notice the singularity at the Heisenberg time r = 1 and the fact that 
cr(oo) exists. 
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However if we smooth out in time R(t), taking for example 

R = — - — P R(r)dr (34) 
Ti - r Jto 

then we get a s elf t- averaging quantity 

JSa < 35 > 

with 

<r = r dra{r) (36) 

Ti - T Jr 

Some numerical work on chaotic billiards [5], in the large window limit, confirm 
this exponential distribution. Integrable billiards show a very different behaviour 
[5]. 

Finally, we would like to mention the fact that Wigner's energy level statistics 
can be obtained for models, where eigenvalues and eigenvectors are correlated. We 
think therefore that the study of quantum dynamics could discriminate between 
such models and those we have considered where they are uncorrelated. 
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